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ABSTRACT 

Apollonian circle packings arise by repeatedly filling the interstices between mutually tan- 
gent circles with further tangent circles. In Euclidean space it is possible for every circle in 
such a packing to have integer radius of curvature, and we call such a packing an integral Apol- 
lonian circle packing. There are infinitely many different integral packings; these were studied 
in the paper 8^. Integral circle packings also exist in spherical and hyperbolic space, provided 
a suitable definition of curvature is used (see ^Oj ) and again there are an infinite number of dif- 
ferent integral packings. This paper studies number-theoretic properties of such packings. This 
amounts to studying the orbits of a particular subgroup A of the group of integral automorphs 
of the indefinite quaternary quadratic form Qt>{w-, x, y, z) = 2{'uP' ^-y"^ + z^) — {w+x+y+ . 
This subgroup, called the Apollonian group, acts on integer solutions QT>iw,x,y, z) = k. This 
paper gives a reduction theory for orbits of A acting on integer solutions to Qviw, x, y,z) = k 
valid for all integer k. It also classifies orbits for all A; = (mod 4) in terms of an extra param- 
eter n and an auxiliary class group (depending on n and k), and studies congruence conditions 
on integers in a given orbit. 
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1. Introduction 



A Descartes configuration is a set of four mutually touching circles with distinct tangents. A 
(Euclidean) Apollonian circle packing in the plane is constructed starting from a Descartes 
configuration by recursively adding circles tangent to three of the circles already constructed 
in the packing; see ^01 ; |Ej ; |Ej i and Figure Q 

Apollonian packings can be completely described by the set of Descartes configurations 
they contain. The curvatures (a, 6, c, d) of the circles in a Descartes configuration satisfy an 
algebraic relation called the Descartes circle theorem. Using the Descartes quadratic form 



Qv{a, b, c, d) := 2{a^ + + + d^) - {a + b + c + df , 



(1.1) 



Figure 1: The Euclidean Apollonian circle packing (—1,2,2, 



^3) 



the Descartes circle theorem states that 



QT){a,b,c,d) = 0. 



(1.2) 



Descartes originally stated this relation in another algebraic form; see |1U| . In order to have 
this formula work in all cases, the curvatures of the circles in a Descartes configuration must 
be given proper signs, as specified in ^U] and 

The form Qx> is an integral quadratic form in four variables, and the equation (|1.2j) has many 
integer solutions. If the initial Descartes configuration generating a (Euclidean) Apollonian 
packing has all four circles with integral curvatures, then it turns out that all the circles in the 
packing have integer curvatures. This integer structure was studied in the first paper in this 
series 8 using the action of a particular discrete subgroup of automorphisms of the Descartes 
quadratic form, which was there termed the Apollonian group A. It is a group of integer 
matrices, with generators explicitly given in §3. The Descartes configurations in a (Euclidean) 
Apollonian packing are described by an orbit of the Apollonian group. That is, the Apollonian 
group acts on the column vector of curvatures of the initial Descartes configuration, and 
generates the curvature vectors of all Descartes configurations in the packing. This fact explains 
the preservation of integrality of the curvatures. The Apollonian group and its relevance for 
Apollonian packings was already noted in 1992 by Soderberg |19j . 

The first paper in this series jSj studied Diophantine properties of the curvatures in Eu- 
clidean integer Apollonian circle packings, based on the Apollonian group action. In ^01 it was 
observed that there exist notions of integer Apollonian packings in spherical and hyperbolic 
geometry. These are based on analogues of the Descartes equation valid in these geometries, as 
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follows. In hyperbolic space, if the "curvature" of a circle with hyperbolic radius r is taken to 
be coth(r), then the modified Descartes equation for the "curvatures" of circles in a Descartes 
configuration is 

Qv{a,b,c,d) =A. (1.3) 

Similarly, in spherical space, if the "curvature" of a circle with spherical radius r is defined to 
be cot(r), then the modified Descartes equation is 

QT,{a,b,c,d) = -A. (1.4) 

There exist integral Apollonian packings in both these geometries, in the sense that the (spher- 
ical or hyperbolic) "curvatures" of all circles in the packing are integral; see the examples in 
Figures |21 and 01 below. Again, the vectors of "curvatures" of all Descartes configurations in 
such a packing form an orbit under the action of the Apollonian group. 

The purpose of this paper is to study Diophantine properties of the orbits of the Apollonian 
group A acting on the integer solutions to the Diophantine equation 

Qv{a,b,c,d) = k, (1.5) 

where is a fixed integer. Integer solutions exist for = or 3 (mod 4), and our results 
on reduction theory in §3 are proved in this generality. Our other results are established 
under the more restrictive condition k = (mod 4), where we will always write k = Am. 
These results cover the motivating cases k = ±4 corresponding to spherical and hyperbolic 
Apollonian packings, as well as the Euclidean case A; = 0. We note that the Apollonian group 
^ is a subgroup of infinite index in the group Aut((5x'7^) of all integer automorphs of the 
Descartes quadratic form Qj). As a consequence, the integral solutions of QT>{a,b,c,d) = k 
(for /c = or 3 ( mod 4)) fall in an infinite number of orbits under the action of A. In contrast, 
these same integral points fall in a finite number of orbits under the full Aut(Qx>, Z)-action. 

In ^ we study the distribution of integer solutions to the equation Qx){a,b,c,d) = 4m, 
determining asymptotics for the number of such solutions to these equations of size below a 
given bound. These asymptotics are obtained via a bijective correspondence between integral 
representations of integers 4m by the Descartes form and representations of integers 2m by 
the Lorentzian form Qc{W,X,Y, Z) = -W"^ + + + Z^, and then applying results of 
Ratcliffe and Tschantz fl7) on integer solutions to the Lorentzian form. 

In ^ we describe the Apollonian group, and apply it to give a reduction algorithm for 
quadruples v = {a,b,c,d) satisfying Qx>(v) = /c in a fixed orbit of the Apollonian group, 
starting from an initial integer quadruple satisfying Qx>(v) = k. This reduction algorithm 
differs in appearance from that given in part I J51 Sect. 3], which treated the case k = 0; however 
we show its steps coincide with those of the earlier algorithm when k = 0. A quadruple obtained 
at the end of this algorithm is called a reduced quadruple. We classify reduced quadruples as 
two types, root quadruples and exceptional quadruples. We show that for each orbit containing 
a root quadruple, the root quadruple is the unique reduced quadruple in that orbit. For each 
fixed k we show that there are at most a finite number of orbits not containing a root quadruple. 
Each of these orbits contain at least one and at most finitely many reduced quadruples, each of 
which is by definition an exceptional quadruple. We show that such orbits can only occur when 
A; > 0, and we call them exceptional orbits. Their existence is an interesting new phenomenon 
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Figure 2: The spherical Apollonian circle packing (0,1,1,2) 




uncovered in this work. They do occur in the hyperboUc Apollonian packing case /c = 4, 
where we show there are exactly two exceptional orbits, labeled by the hyperbolic Descartes 
quadruples (0, 0, 0, 2) and (—1, 0, 0, 1), respectively. We give a geometric interpretation of these 
exceptional hyperbolic Apollonian packings. 

As examples, Figure|21shows the spherical {k = —4) Apollonian packing with root quadruple 
(0,1,1,2), after stereographic projection. Figure shows the hyperbolic {k = 4) Apollonian 
packing with root quadruple (—1,1,1,1), using the unit disc model of the hyperbolic plane. 
For more on the geometry of Apollonian packings in spherical and hyperbolic space consult 

m- 

In 311 we count the number A''i.oot(4m, — n) of root quadruples a = (a, b, c, d) with Qv{^) = 
4m having smallest member a = —n. We give a formula for this number in terms of the 
class number of (not necessarily primitive) integral binary quadratic forms of discriminant 
— 4(n^ — m) under GL(2, Z)-equivalence. The class number interpretation applies only when 
n? > m. For certain m > there may exist root quadruples with parameters < m, 
and information on them is given in Theorem 13.31 in ^ For the hyperbolic case the only 
value not covered is n = — 1. Associated to it is an infinite family of distinct root quadruples 
{(— 1, 1, c, c)| c > 1}. FigureOlpictures the Apollonian packing with root quadruple (—1, 1, 1, 1). 

The class number equivalence allows us to derive good upper bounds for the number of 
such root quadruples with smallest value a = —n as n — > oo, namely 

A'root(4m, -n) = O (n(logn)(loglogn)^) . 

If m = ±1, then these bounds are interpretable as bounds on the integral spherical and 
hyperbolic packings with boundary an enclosing circle of (spherical or hyperbolic) curvature 
n. 

In ^ we investigate congruence restrictions on the integer curvatures which occur in a 
packing, for the case of spherical and hyperbolic packings, i.e., k = ±4. We show that there 
are always non-trivial congruence conditions modulo 12 on the set of integers that can occur 
in a given packing. It seems reasonable to expect that such congruence conditions can only 
involve powers of the primes 2 and 3 and that all sufficiently large integers which are not 
excluded by these congruence conditions actually occur. 

Various open problems are raised at the end of ^21 33 In addition, in the case of 

Euclidean Apollonian packings, it is known that there are Apollonian packings which have 
stronger integer properties, involving the centers of the circles as well as the curvatures. These 
were studied in Graham et al. 0, [HI, and j7j. It remains to be seen if there are analogues of 
such properties (for circle centers) in the spherical and hyperbolic cases. 

Acknowledgments. Much of this work was done while the authors were at AT&T Labs- 
Research, whom the authors thank for support. N. Eriksson was also supported by an NDSEG 
fellowship. The authors thank the reviewer for helpful comments. 

2. Integral Descartes Quadruples 

The Descartes quadratic form Qx> is the quaternary quadratic form 

Qviw, X, y, z) = 2{n? + + + z^) - {w + x + y + z f. (2.1) 
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Its matrix representation, for v 



IT 



[w,x,y,z\- , IS 
/I 



V- 



1 



-1\ 

-1 
-1 

1 / 



V. 



This form is indefinite with signature (+, +, +, — ). We consider integral representations of an 
integer k by this quaternary quadratic form. For an integer quadruple v = {w,x,y,z) G Z'*, 
Qx>{w,x,y, z) = or 3 (mod 4), according to the parity oi w + x + y + z. We define the 
Euclidean height H{v) of a (real) quadruple v = (w, x, y, z) to be 

H{^,) := (u;2+x2+y2 + ^2)l/2_ 

Notice that we reserve the usual notation, |?;|, for another meaning in ^ 

We relate integer representations of even k of the Descartes form to integer representations 

2 



of t of the Lorentz form 



QciW, X, Y, Z) = -W'^ + X'^ + Y^ + Z^. 

in Lemma l2 . 1 1 b elow . The relation only works for even k, so in this section we consider only 
the case /c = (mod 4), so that 

Qviw,x,y,z) = k = im. (2.2) 
The Lorentz form has the matrix representation, for v = [W,X,Y, Z]'^ , 



/-I 




Vo 





1 
1 







V 



The Descartes form and Lorentz form are related by 



2Qc = JoQv^o , 



(2.3) 



where 



/ 1 
1 

1 

V 1 



1\ 

-1 
-1 

1 / 



Notice that Jg = I. 

Let Nc{k,T) count the number of integer representations of k by the Lorentz form of Eu- 
clidean height at most T and let Nx>{k, T) denote the number of solutions of (|2.2|) of Euclidean 
height at most T. 

Lemma 2.1. The mapping (W, X,Y, Z)'^ = 3Q{w,x,y,z)'^ gives a bijection between real so- 
lutions of Qx){w, x,y, z) = 4r and real solutions of Qc{W, X,Y, Z) = 2r for each r G M, and 
this bijection preserves the Euclidean height of solutions. In the case where r = m £ Z this 
bijection restricts to a bijection between integer representations of Am by the Descartes form 
and integer representations of 2m by the Lorentzian form. In particular, for all integers m, 

NT,iAm,T) = Nc{2m,T). 
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Proof. There is a bijection on the level of real numbers, by (|2.3() . since Jq is invertible, and 
Qviw, X, y, z) = 4r if and only if Qc{W, X, Y, Z) = 2r. A calculation using Jo = I shows that 

H{w, X, y, z) = + x^ + y^ + x^ = + + + = H{W, X, Y, Z). 

Now suppose that r = m £ Z. The mapping takes integral solutions {w, x, y, z) to 
integral solutions {W, X,Y, Z) because all integral solutions of Qx>{w, x,y, z) = Am satisfy 
w + x + y + z = (mod 2), as can be seen by reducing the Descartes equation (|2.2I) modulo 
2. In the reverse direction integral solutions go to integral solutions because all solutions of 
Qc{w, X, y, z) = 2m have W + X + Y + Z = Q (mod 2), similarly. □ 

Representation of integers by the Lorentz form has been much studied; see Radcliffe and 
Tschantz jl7j . Their results immediately yield the following result. 

Theorem 2.2. For each nonzero integer m, there is an explicitly computable rational number 
c{2m) such that 

Nv{Am,T) = c{2m) + o{T^), (2.4) 

in which 

n=0 ^ ' 

We have c(-2) = | and c(2) = |. 

Proof. From Lemma l2. li the problem of counting integral solutions to Qti{w, x,y,z) = Am of 
height at most T is equivalent to counting the number of integral solutions to (^^(VF, X, Y, Z) = 
2m of Euclidean height at most T. We use asymptotic formulae of Radcliffe and Tschantz for 
the number r{n, k, T) of solutions to the Lorentz form —Xq+X^ + • • ■ + X'^ = k with |Xo| < T, 
for any n > 2 for any fixed nonzero k. We take n = 3, and their formula for A; < in |171 
Theorem 3, p. 505] is 

r(3, k, T) = M:^5(3, k)T^ + 0(t3/2), ^2.5) 

in which Vo^S*^) = vr and (5(3, A;) is a density constant for the representation in the sense of 
Siegel USl- We need asymptotics for the number of solutions s{n,k,T) to the Lorentz form 
QciW,X,Y, Z) = k of Euclidean height below T, which is given by 



s{n,k,T) :=r ( n,k, ^1 ^{T^ - k) 



The asymptotic formula (|2.5j) yields 



s{3, k, T) = ^ (5(3, 2m)T^ + 0{T^/^). (2.6) 
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Their result ^1 Theorem 12, p. 518] evaluates the constant (5(3, k) in terms of the value of an 
L-function, as 



in which each Sp{3, k) is an effectively computable local density which is a rational number. 
We now set k = 2m, and obtain a formula of the desired shape (|2.4|) . with 



c(2m) 



1 / 6p{3,2m) 

n(i-(^)p- 



Radcliffe and Tschantz |17[ Theorem 4, p. 510] also obtained an asymptotic formula for 
solutions to the Lorentz form with A; > 0, which gives 

r{3,k,T) = ^^^5{3,k)T^ + o{T^). 

We derive a similar asymptotic formula for s{n,k,T) (without explicit error term), by setting 
k = 2m, and proceeding as above. 

The particular values c(— 2) = | and c(2) = | are derived from Table II in Radcliffe and 
Tschantz p. 521], using k = ±2. □ 



Remarks. (1) The case A; = exhibits a significant difference from the pattern of Theo- 
rem In ^ Theorem 2.1] it was shown that 



7r2 



4L(2,x-4) 

That is, the coefficient c(0) = f is transcendental. 

(2) Theorem 12.21 implies that the number of hyperbolic Descartes quadruples of Euclidean 
height below T is asymptotically | the number of spherical Descartes quadruples of height 
below r, as r — > 00. That is, 

iV.(4.T) ^ 5 

T^oo Nv{-4.,T) 3 ^ ^ 



3. Reduction Theory and Root Quadruples 

The results of this section apply to the general case 

Qv{w,x,y,z) = k, (3.1) 

with A: = or 3 (mod 4) an integer. We introduce the Apollonian group, and present a reduc- 
tion procedure which, given an integral Descartes quadruple, uses the action of the Apollonian 
group to transform it to a reduced quadruple, which is minimal according to a certain measure. 
The reduced quadruples then classify integral Apollonian circle packings, generally uniquely, 
but up to a finite ambiguity in some exceptional cases. 
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3.1. Apollonian Group 



In |H] and ^Hl Apollonian circle packings were specified in terms of the set of Descartes con- 
figurations they contain. This set of configurations completely describes the packing, and it 
was observed there that, in a suitable coordinate system, they form as the orbit of any single 
Descartes configuration in the packing under the motion of a discrete group, the Apollonian 
group. 

Definition 3.1. The Apollonian group A = (Si, S2, S3, S4) is the subgroup of GL(4, Z) gen- 
erated by the four integer 4x4 matrices 



Si = 



/-I 


2 


2 


2^ 




A 








o\ 





1 










2 


-1 


2 


2 








1 





, S2 = 








1 

















Vo 








1/ 


/I 








^\ 




/I 








o\ 





1 













1 








2 


2 


-1 


2 


, S4 = 








1 





\o 








V 




\2 


2 


2 


-1/ 



The Apollonian group is a subgroup of the integer automorphism group Aut((5D, Z) of the 
Descartes form under congruence, which is the largest subgroup of GL(4, Z) which leaves Qd 
invariant; that is, 

U'^QpU = Qv, for all U G ^, (3.2) 

a relation which needs only be checked on the four generators Sj. In particular it preserves 
the level sets (5o(v) = k. It acts on all (real) Descartes quadruples v = (a, b, c, d)^ (viewed as 
column vectors) via matrix multiplication, sending v to Uv. 

The group Ant{QD,Z) is very large, and can be put in one-to-one correspondence with 
a finite index subgroup of the integer Lorentzian group 0(3, 1,Z) (using the correspondence 
(|2.3j) in §2). The Apollonian group A is of infinite index in Aut^Qxi,'^), and it has infinitely 
many distinct integral orbits for each integer A; = or 3 (mod 4), the values for which integral 
solutions exist. 

Geometrically, the generators Sj in A correspond to inversions with respect to a circle 
passing through three intersection points in a Descartes configuration. This inversion gives 
a new Descartes configuration in the same packing, leaving three of the circles fixed; see 
Section 2]. In particular the group can be viewed as acting on column vectors v = [a, 6, c, d]^ 
giving the "curvatures" in a Descartes configuration, converting it to the curvatures of another 
Descartes configuration in the packing. The generators Sj allow us to move around inside the 
packing, from one Descartes configuration to a neighboring configuration. From (|2.1j) and (|3.1|) 
we see that if the curvatures a, b, c of three touching circles are given, then the curvatures of 
the two possible circles which are tangent to these three satisfy 

d,d' = a + b + c± y4(a6T"krracyTA;. 
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Therefore d + d' = 2{a + b + c). If we begin with circles with curvatures a, b, c, d, the curvature 
d'of the other circle which touches the first three is 

d' = 2{a + b + c)- d. 

The Apollonian group action corresponding to this sends (a, b, c, d) to (a, b, c, d'), which is the 
action of S4. In particular, if we start with an integer quadruple, this procedure will only 
create other integer quadruples. This group is studied in Aharonov Graham et al. [3], [S] 
and Soderberg [T^ . 

Definition 3.2. An integral Descartes ensemble ^[v] is an orbit of the Apollonian group A 
acting on an integral quadruple v = {w,x,y, z)'^ G Z^. That is, Al'v] := {Uv : U E A}. 

All elements x G .4[vo] are integer quadruples satisfying 

(3d(x) = (5d(vo) = k. 

We are most interested in the case L{v) :=w + x + y + z = ( mod 2), where k = 4m. For the 
cases m = —1,1 such ensembles give the set of Descartes quadruples in an integral spherical 
(resp. hyperbolic) Apollonian circle packing. 

This section addresses the problem of classifying the orbits of the Apollonian group acting 
on the set of all integral solutions Qvi'v) = k for a fixed value of k. 

3.2. General Reduction Algorithm 

We present a reduction procedure which, given an integral quadruple vq, finds an element of 
(locally) "minimal" size in the orbit ^[vq]. We will show that every orbit contains at least one 
and at most finitely many reduced elements. In most cases the "minimal" quadruple is unique 
and is independent of the starting point vq of the reduction procedure, and this is the case 
for root quadruples defined below. However when k > there sometimes exist exceptional 
integral Descartes ensembles containing more than one reduced element, and the outcome of 
the reduction procedure depends on its starting point in the orbit. We show that for each fixed 
k there are in total at most finitely many such exceptional orbits. 

The general reduction procedure given below greedily attempts to reduce the size of the 
elements in a Descartes quadruple v := {a,b,c,d)'^ by applying the generators Sj to decrease 
the quantity 

|v| := \a\ + |6| + |c| + 

General reduction algorithm. 

Input: An integer quadruple (a, b, c, d). 

(1) Order the quadruple so that a < b < c < d. Then test in order 1 < i < 4 whether some 
Si decreases |v| := \a\ + \b\ + |c| + \d\. For the first Sj that does, apply it to produce a new 
quadruple, and continue. 

(2) If no Sj strictly decreases |v|, halt. Declare the result a reduced quadruple. 



10 



This reduction procedure is slightly different from the reduction algorithm used in part 
I IS] for the case k = 0. The part I reduction procedure applied only to quadruples with 
L(a) = a + b + c + d > 0, and tried to decrease the invariant L{a.) = a + b + c + d at each step, 
halting if this could not be done. Although the general reduction algorithm uses a different 
reduction rule, one can show in this case that it takes the identical sequence of steps as the 
reduction algorithm in part I; see the remark at the end of §3.3. Thus the general reduction 
algorithm can be viewed as a strict generalization of the algorithm of |H] . 

We will classify reduced quadruples as either root quadruples or exceptional quadruples, as 
defined in the following theorem. 

Theorem 3.1. The general reduction algorithm starting from any nonzero integer quadruple 
a = (a, b, c, d) always halts in a finite number of steps. Let the reduced quadruple at termination 
be ordered as a = (o, b, c, d) with a < b < c < d. Suppose L{a) = a + b + c + d>0. Then 
exactly one of the following holds: 
(i) The quadruple a satisfies 

a + b + c>d>0. (3.3) 

In this case we call a a root quadruple, 
(a) The quadruple a satisfies 

a + b + c<0<d. (3.4) 

In this case we call a an exceptional quadruple. 

If a is a reduced quadruple such that L{a) < 0, then a* := (— d, — c, —6, — o) is a reduced 
quadruple with L{a*) > 0. We call such a a root quadruple (resp. exceptional quadruple) if 
and only if a* is a root quadruple (resp. exceptional quadruple). 

Proof. The integer-valued invariant |a| := \a\ + \b\ + |c| + \d\ strictly decreases at each step of 
the reduction algorithm. It is nonnegative, so the process stops after finitely many iterations. 

Now suppose that a is reduced and L{a) = a + b + c + d > 0. We must show that (i) or (ii) 
holds. Now d > because a is nonzero. The condition that a be reduced requires |S4a| > |a|, 
and since S4a = {a,b,c,d') with d' = 2{a + b + c) — d, this condition is \d'\ > \d\ = d. Now 
d' > d gives a + b + c'> d which is (i), while d' < —d gives a + 6 + c < which is (ii). □ 

The definition of root quadruple formulated in Theorem 13.11 is not identical in form to the 
definition of root quadruple used in part I [H] for the case k = 4m = 0. However they are 
equivalent definitions, since the reduction algorithms take identical steps, as remarked above. 

3.3. Root Quadruples 

We first show that the reduction algorithm is well behaved for all orbits ^[v] containing a root 
quadruple. 

Theorem 3.2. Let A[v] be a nonzero integer orbit of the Apollonian group, and suppose that 
A[v] contains a root quadruple a. Then: 
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(1) The quadruple a is the unique reduced quadruple in A\v\. 

(2) All values L{x) for x £ A[v] are nonzero and have the same sign. If this sign is positive, 
then the general reduction algorithm starting from any x € ^[v] applies only S4 until the root 
quadruple is reached. 

Proof. It suffices to prove the result in the case when the root quadruple has L{a) = 
a + b + c + d > 0, because the case L(a) < can be treated using a*. Without loss of generality 
we may reorder the quadruple a < b < c < d (since A is invariant under conjugation by 
elements of the symmetric group on 4 letters, represented as permutation matrices). Now the 
root quadruple condition ()3.3|) yields 

L(a) = a + 6 + c + d > a + 6 + c > d > 0. (3.5) 

This implies that 

a + b>d-c>0, (3.6) 
where the last inequality follows from the ordering. 

Claim. Let b G A[v] be a Descartes quadruple and let Po-b = {w, x, y, z) be its coordinates 
permuted into increasing order w < x < y < z. Then these coordinates satisfy: 

Property (PI): w + x > 0. 

Ifh does not equal the given root quadruple a then these coordinates satisfy: 

Property (P2): ^<w + x + y<z. 



To prove the claim, we observe that every quadruple b 7^ a in A\v] can be written as 
b = Sj^Sj^_^ • • • Sj^a with each Si- G {Si, S2, S3, S4}, for some /c > 1, and with Sj^ 7^ 
2 < j < fc, since all S| = /. We proceed by induction on the number of multiplications k. 

The case /c = is the root quadruple a, which satisfies property (PI) by 1)3. For the base 
case A; = 1 we must verify that properties (PI) and (P2) hold for S^a for 1 < j < 4, where a 
is the given root quadruple. Without loss of generality we may assume a = (a, 6, c, d) has the 
ordering a < b < c < d] a, permutation of a merely permutes the four values S^a. Consider 
first Sia = (a', 6, c, d). We assert that the ordering of the elements of Sia is 

b<c<d<a'. (3.7) 

To establish this we first show 

a' ■.= 2{b + c + d)-a = {b + c + d) + {b + c + d-a)>b + c + d. (3.8) 

This inequality holds since 6 — a > from the ordering while d > and c > from (|3.5)) . This 
yields a'>b + c + d> {a + b) + d > d, which establishes (|3.7() . Property (PI) for Sia asserts 
that 6 + c > 0, and this holds by (|3.6j) . Property (P2) for Sia asserts that Q<b + c + d<a\ 
and the left inequality follows from b + c + d > {a + b) + d > 0, using both 1)3. 5|) and 1)3. 6() . 
while the right inequality follows from ()3.8j) . This proves the base case for Sia. The proofs 



12 



of the base cases S^a for j = 2,3 are essentially the same as for Sia. In these cases the one 
curvature that changes becomes strictly larger than d. Finally we treat = (a, b, c, d'). The 
ordering that holds here is 

a<b<c<d'. 

Indeed we have 

d' = 2{a + b + c)-d={a + b + c) + {a + b + c-d)>a + b + c>d, (3.9) 

by H3.5() . If d = d' we stay at the root quadruple, so we may assume d' > d in what follows. 
Property (PI) for S4a asserts a + b > 0, and this holds by (|3.6|) . Property (P2) for S4a asserts 
0<a + b + c<d'. Now if a + 6 + c = d then d' = d, so assuming d' > d we must have 
a + b + c > d. Using this fact, we see that the inequalities in (|3.9|) are both strict, which gives 
the right inequality a + b + c < d'. The left inequality a + 6 + c > holds by ()3.5|) . Thus 
properties (PI) and (P2) hold for S4a as well, or else we stay at the root quadruple. This 
verifies (PI) and (P2) in all cases, and shows that if b = S^a ^ a, then the largest coordinate 
in b is unique and is the coordinate in the j-th position. 

For the induction step, for A; + 1 given k, we suppose that b = Sij.Sjj._-^ • • • Sj^a with b 7^ a, 
k > 1, and we include in the induction hypothesis the assertion that the largest coordinate of 
b is unique and occurs in the ifc-th coordinate position. Let b' = {w, x, y, z) = P^h denote its 
rearrangement in increasing order w < x < y < z, where Po- is a 4 x 4 permutation matrix. The 
induction hypothesis asserts that (PI) and (P2) hold for b', together with the strict inequality 
y < z. The induction step requires showing that properties (PI) and (P2) hold for Sjb' for 
1 < j < 3, corresponding to a product of length k + 1, and that the new coordinate formed in 
Sjh' is its strictly largest coordinate. The case of multiplication by S4 is excluded because it 
corresponds to choosing Sj^.^-^ = Sj^. . Indeed S4 changes the largest coordinate, which Sj^. did 
at the previous step by the induction hypothesis, so that Si^^-^ := P~^S4Po- = Sj^.. The first 
case is Sib' = {w' ^x^y^ z). We assert that 

x<y<z<w'. (3.10) 

Indeed x, y, z > by property (PI) for b' and 

w' = 2{x + y + z)— w = {x + y + z) + {x + y + z — w)>x + y + z, 

since property (P2) for b' gives x + y + z>w + 2x + 2y>w. Next property (PI) for b' gives 

X + y + z > {w + x) + z > z > y > 

whence w'>x + y + z>z and (|3.10() follows. Property (PI) for Sib' asserts x + y > 0, 
which holds since x + y>w + x>{)^ using Property (PI) for b'. Property (P2) for Sib' 
asserts 0<x + y + z< w' , both inequalities of which were shown above. The second case is 
S2b' = {w,x',y,z). We assert that 

w <y < z <x', (3.11) 

where we must show x' > z. Now x' = 2{w + y + z) — x = {w + y + z) + {w + y + z — x) > w + y + z 
since property (P2) for b' gives w + y + z > 2w + x + 2y > x. Next w + y>w + x>Q holds 
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by property (PI) oih',sow + y + z>z>0 whence x' > z, and follows. Property (PI) 

for S2b' now asserts w + y > 0, which holds since w + y + z > z. Property (P2) for S2b' now 
asserts < w + y + z < x' , and both of these inequalities were verified above. Finally, the 
proof for Ssb' = {w,x,y',z) is similar, where one shows w < x < z < y'. This completes the 
induction step, and the claim follows. 

To continue the proof, we observe that for any non-root quadruple b' property (P2) implies 
that L{h) = {w + X + y) + z > 0. Thus all quadruples in the packing have L{h) > 0. We 
next show b is not reduced. Indeed property (P2) gives S4b = {w, x, y, z') with < z' = 
2{w + x + y) — z<z so that |S4b| < |b|. We conclude that the root quadruple is the unique 
reduced quadruple, which proves assertion (1). 

To verify assertion (2), we must show that |Sjb| > |b| for 1 < j < 3, while |S4b| < |b|. 
The last inequality is already done. Now S^h = {w,x,y',z) with y' = 2{w + x + z) — y > 
2{w + x) + 2{w + x + y)— y>y has jSsbl > |b|. The argument showing |Sjb| > |b| if j = 1,2 
is similar. This proves (2). □ 

For Qxi{ai) = k with k fixed there are infinitely many distinct root quadruples. The next 
result determines various of their properties. 

Theorem 3.3. (1) Any root quadruple a = (a, b, c, d) with Qx>{sl) = k and a < b < c < d and 
L(a) > satisfies 

0<b<c<d. (3.12) 

In addition, if k > 0, then 

a < 0, (3.13) 

and if k < then 

a < (3.14) 

(2) For each pair of integers {k,—n) there are only finitely many root quadruples a with 
(5©(a) = k and a = —n, with the exception of those pairs (k, —n) = {AP, —I) with I > 0. In 
each latter case there is an infinite family of root quadruples {{—I, l,c,c): c > max(/, 1)}. 

Proof. (1) The condition 1)3. 3() implies a-\-b>d — c>0, so that b > 0, and the ordering 
gives pT2ll . 

We view the Descartes equation Qx>{^) = as a quadratic equation in a, and solving it 
gives 

a = b + c + d± a/ 4{bc + bd + cd) + k 

The ordering and 1)3. 12() give a < b < b + c + d, which shows that the minus sign must be taken 
in the square root, and 

a = b + c + d- V4(6c + bd + cd) + k (3.15) 
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Using H3.3() we have 

2{bc + hd + cd) = {be + d{b + c)) + {be + bd + ed) 
> {b^ + d{d - a) + {be + bd + c^) 
= b"^ + + d^ + [be + d{b - a)], 

which gives 

4(6c + bd + ed)>{b + e + df + [be + d{b - a)]. 
As a consequence (|3.15() yields 



a<{b + e + d)- {b + e + df + k. (3.16) 
For A: > this immediately gives 



a<b + e + d- ^{b + c + dY + 1 < 0, 

which is (TTm . 

For A; < we must have w := b + c + d> \/\k\ so that the square root above is real, whence 
(|3.16|) gives 



a < max [w — yuP^+k] = 

w>^J\k\ 



which is (|3.14j) . Indeed the function f{w) = w — Vw"^ + k on this domain is maximized at its 
left endpoint w = \/\k\, since 

f'{w) = 1 - -j=^== < for ^ > 

Root quadruples with a > do occur for some negative k. For example v = (2, 3, 4, 7) is a 
root quadruple for k = Qt>{^) = —96. 

(2) We take the values Qx'(a) = k and a = —n as fixed. To show finiteness, it suffices to 
bound c above, in terms of k and n. The conditions a + b + e > d and Q < b < c < d yield 
< d < 3c, < 6 < c and -2c < a < c. 

We have d = a + b + e — x for some x > 0, and substituting this in ()3.15j) gives 



a = a + 26 + 2c-x - ■v/4(6 + c)(a + 6 + c - x) + 46c + A:. 
This simplifies to 

= 26 + 2c + X - \/(26 + 2c - xY + (-x^ + 4a6 + 4ac + 46c + A;), 
so we must have 

x'^ = Aab + Aae + Abc + k. (3.17) 

Now d>e+{a + b — x) > c gives 

a + 6 > X > 0. 
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Squaring the left inequality and rearranging gives + 2ab + 6^ > = 4ab + 4ac + 46c + k, 
which yields 

{a-bf >4:{a + b)c + k (3.18) 

We have c > 6 > 0, and holding a and k fixed and letting b grow, one sees that the right side 
grows at least as fast as 46^ while the left side grows like 6^, so one concludes that b is bounded 
above. One can now check that ()3.18|) implies 

b < 3\a\ + 

This bound shows that the left side of (|3.18|) is bounded above, yielding 

4(a + 6)c< (4|a|)2 + |A;|. 

We have a + 6 > x > 0, and whenever a + 6 > then this inequality bounds c above, with 
c < + \k\, as was required. 

There remains the case a+b = x = 0. Letting a = —I, we have b = I and d = a+b+c—x = c. 
Thus we have the quadruple a = (—l,l,c,c) with the requirement c> I > 0. All of these are 
root quadruples except for / = c = which gives the excluded value (0, 0, 0, 0). Here (5(a) = 'iP, 
and for each / > we obtain the given infinite family of root quadruples. □ 

Remark. For k > the upper bound a < \/\k\ in Theorem 13. 3r i) is not sharp in general. For 
the case k = —4 of spherical root quadruples, this upper bound gives a < 2, but all spherical 
root quadruples satisfy the stronger bound a < 0. To see this, observe that a = 2 can only 
hold with b + c + d = 2in (|3.16jl , and since b < c < d are nonnegative integers we must have 
6 < so a < 0, a contradiction. Similarly a = 1 gives l<6 + c + d<fin (nOT)!) . with the 
same contradiction. The case A; = does occur with the spherical root quadruple (0,1,1,2) 
pictured in Figured 

Remark. We conclude this subsection by showing that the notions of reduction algorithm 
and root quadruple used here agree for k = with those used in part I The part I 
reduction procedure applied only to quadruples with -L(a) = a + b + c + d > 0, and tried to 
greedily decrease the invariant L{sl) = a + b + c + d at each step, halting if this could not 
be done. Although the general reduction algorithm uses a different reduction rule, in this 
special case it takes the identical series of steps as the reduction algorithm in part I, so defines 
the identical notion of root quadruple. Indeed, suppose k = and L{a) > 0. We will need 
a result from Theorem I3.4f 3) below, which shows that there are no exceptional quadruples 
for k = 0. so the general reduction algorithm always halts at a root quadruple. Comparing 
Theorem I3.2r 2) with Lemma 3.1 (iii) of ,5, shows all the steps are identical. Namely, after 
reordering the curvatures in increasing order, both algorithms always apply S4. Finally we 
check that the halting rules coincide, i.e., that |a| is minimized where L{aL) is minimized. For 
both algorithms the minimal element in a root quadruple is non-positive, with the other three 
elements nonnegative (see Theorem 13. 3r i) and |j5i Lemma 3.1(3)]). Furthermore the first non- 
positive element encountered in either algorithm is never changed subsequently. Since the two 
algorithms coincide up to this point, and afterwards the two invariants remain in the fixed 
relation |a| = L{a.) + 2\a\, for fixed |a|, their halting criteria coincide. 
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3.4. Exceptional Quadruples 



We next consider orbits ^[v] that contain exceptional quadruples. 

Theorem 3.4. Let A[v] be a nonzero integer orbit of the Apollonian group. 

(1) If A\v] contains an exceptional quadruple, then it contains elements vi, V2 with i(vi) > 
and -L(v2) < 0. 

(2) There can be more than one exceptional quadruple in A\v]. 

(3) Exceptional quadruples can exist only for k >1. For each A; > 1 there are only finitely 
many exceptional quadruples with Qt>{^) = k, clU with H^a)'^ < 2k'^ . 

Proof. (1) Let a = {a,b,c,d) be an exceptional quadruple in ^[v], with a < b < c < d. It 
suffices to treat the case L(a) = a + b + c + d>0, since we can otherwise apply the argument 
to a*. 

Now Theorem 13.11 gives a + 6 + c < and d > 0. Take vi = Sia = (a', b, c, d)-^ and 
L(vi) = 3(a + b + c + d) — 4a > if a < 0, while if a = then a = b = c = and so d > and 
L(vi) > in this case. Next, take V2 = S4a = (a, b, c, d'). Then L(v2) = 3(a + b + c) — d<0 
if d > 0, while if d = then a = b = c = d = 0, a contradiction. 

(2) The quadruple a = (—2, —3, —5, 12) with |a| = 22 and Qvi^) = k = 360 is an ex- 
ceptional quadruple. Its neighbors are (—2, —3, —5, —32), (—2, —3, 9, 12), (—2, 7, —5, 12), and 
(6, —3, —5, 12). Now (—2, —3, 9, 12) has a neighbor b = (—2, —3, 9, —4) which is also an excep- 
tional quadruple. 

(3) Let k be fixed. Let a = (a, b, c, d) be an exceptional quadruple ordered with a < b < 
c < d. By hypothesis 

|Sj(a)|2 > |a|2 for 1 < j < 4. (3.19) 
We will show that necessarily k > 1 and that 

F(a)2 = 0^ + 52 + c2 + d2 < 2A;2. (3.20) 

The minimality property (|3.19|) and the bound (|3.2()j) are both invariant under reversing all 
signs, so by multiplying by —1 if necessary we may suppose that L(a) = a+b+c+d > without 
loss of generality. Since {a,b,c,d) 7^ (0,0,0,0) we have d > 0. The condition |S4ap > |ap 
requires that either (i) a+b+c > dor (ii) a + b+c < 0. By Theorem 13.11 exceptional quadruples 
correspond to case (ii), and those with L(a) > are characterized by the condition 

d>0>a + b + c. 

This condition necessarily implies a < 0. We consider three exhaustive cases. In each case, 
we change notation so that a, b, c are nonnegative; for example, if a < < b < c < d, then we 
multiply a by —1 so that a = (—a, b, c, d). 

Case 1. -a < -b < -c < < d. 

We have 

k = QT){ai) = + b^ + c^ + d^ -2ab-2ac-2bc + 2ad + 2bd + 2cd 
= + b'^ + + d^ + 2b{d - a) + 2a{d - c) + 2c{d - b) 
> + b^ + c^ + d^ = H{af. 
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Since a / (0, 0, 0, 0) we infer k > 1, and also (|3.2U|) holds in this case. 

Case 2. -a<-b<0<c<d. 
We have 

k = Qv{ai) = a' + + + <f -2ah + 2ac + 2hc + 2ad + 2hd-2cd 
= [a-bf + [c- df + 2ac + 2hc + 2ad + 2bd. 

Every term on the right is nonnegative, which implies A; > 0. This gives the additional 
information (since a> b and d > c) that 

b < a < b + \fk and c < d < c + \Pk. 

as well as ac, &c, ad, bd < ^. If either of a, 6 is nonzero then each of c, d is at most |. Similarly 
if either of c, d are nonzero then a, b are at most |. Now suppose both a = and 6 = 0. Then 
k = (c — d)^, and there is a reduction step that takes (0, 0, c, d) i— > (0,0, c, 2c — d), whence 
|2c — d| > d, so that either c = d, or c = 0. The case c = d gives a root quadruple, which is 
excluded, and in the remaining case (0, 0, 0, d) we have k = d' > 0, and then c,d< ^/k, and in 
all these subcases 

Hia)' =a'^ + b' + + d' Kk' 

holds. The argument if c = d = is similar. In all these cases we infer that k > and that 
(TOni) holds. 

Case 3. -a<0<b<c<d. 

We have 

k = Qv{a) = a' + b' + 0^ + d' + 2ab + 2ac-2bc + 2ad-2bd-2cd 

= {a + bf + {c-df + 2{a-b){c + d) (3.21) 

In the exceptional case we must have —a + b + c < 0, so that a > b + c. This gives a > b, 
whence all terms on the right side of 1)3.21(1 are nonnegative, and one is strictly positive so 
A; > 0. Also (|3.21() gives a,b < Vk, and using a > b + c it also gives 

k>{c-df + 2c(c + d) > 3c' + d'. 

It follows that 

Hia)' = a' + b' + c' + d' <2k < 2k', 

and (pnnil holds. □ 

On comparing Theorem l3.3l and Theorem l3.4l we see that: An integer Apollonian group orbit 
A[v] contains an exceptional quadruple if and only if it contains quadruples with -L(xi) > 
and L{x2) < 0. 

We conclude this section by determining the exceptional quadruples for spherical and hyper- 
bolic Apollonian packings, i.e., for k = ±4. There are no exceptional quadruples for spherical 
packings (k = —4) by Theorem 13.41 (3), so we need only treat the hyperbolic case. 
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Theorem 3.5. For k = A the exceptional orbits A\v\ are exactly those containing an element 
that has a zero coordinate. There are exactly two such orbits, one containing two exceptional 
quadruples (0,0, 0, ±2), and the other containing the exceptional quadruple (—1,0,0,1). 

Proof. By Theorem I3.2r 3). there are finitely many exceptional quadruples for k = A. From 
H3.2U() . we see that if {a,b,c,d) is an exceptional quadruple, then a'^ + b'^ + + d? < 32. A 
computer search shows that there are exactly 7 quadruples with a < b < c < d satisfying this 
bound along with a + b + c + d>0. They are 

{(-1,0, 0,1), (0,0, 0,2), (-1,1, 1,1), (-1,1, 2, 2), (-1,1, 3, 3), (0,0, 1,3), (0,0, 2,4)}. 

It is easy to check that the first five are all reduced while the last two are not. Furthermore, 
the third, fourth and fifth quadruples are root quadruples, so may be discarded. Thus only 
(— 1, 0, 0, 1), (0, 0, 0, 2) remain, and these are exceptional. The orbit of (0, 0, 0, 2) also contains 
the reduced quadruple (—2,0,0,0); see Figure^ The packing with exceptional quadruple 
(—1,0,0, 1) is shown in Figure El 

If fact, if .4[v] does contain a circle of curvature 0, then it necessarily contains an exceptional 
quadruple. Setting a = 0, the orbit contains a integer quadruple (0, b, c, d) with 

F(b, c, d) := b^ + c^ + d^ - 2{bc + bd + cd) = 4. (3.22) 

Without loss of generality, assume that the elements are ordered |6| < |c| < \d\ and that 6 < 0. 
We claim that necessarily 6 = 0. 

Case 1. 6 < 0, c> 0, d > 0. 
Now 

f{b,c,d) = b^ + c^ + d^ + 2\bc\+2\bd\-2\cd\ 
= b^ + (c-df + 2\bc\+2\bd\ 

If \b\ > 1, then |c|, \d\ > 1 also and F{b, c, d) > 1^+0^ + 2 + 2 = 5 which gives us a contradiction. 
So we must have 5 = 0. 

Case 2. 6 < 0, c < 0, d > 0. 
Here |6| < |c| < \d\ gives \bd\ > \bc\, so 

F{b,c,d) = b"^ + c^ + d^ + 2\bd\+2\cd\-2\bc\ 
> b"^ + + d'^ + 2\cd\ 

If 1^1 ^ 1) then so are |c| and \d\, so F{b,c,d) > 5, a contradiction. Thus 6 = 0. 
Case 3. b,c,d< 0. 

First, apply the reduction algorithm to (0, b, c, d). Notice that the first coordinate always stays 
because Si can never make it any smaller. The algorithm halts when S4 can no longer reduce 
\b\ + |c| + \d\, that is, when > \d\. Combining this with d + d' = 2(6 + c), we get \d\ > \b + c\. 
Therefore 

6^ + + < 6^ + + (6 + cf = 26^ + 26c + 20^. 
Since 6c, 6d, cd > 0, we have bc + bd + cd> |6p + |6c| + |cp. Thus, 
F(6, c, d) = 6^ + + d^ - 2(6c + 6d + 6d) < 



19 



and hence cannot satisfy ()3.22|) . a contradiction. So there are no solutions in this case. 

The claim 6 = is proved, so any quadruple with one element zero is contained in the 
same packing as some (0,0, c, d). Now solve 2(c^ + d?) = (c + d)^ + 4 and we see that we 
must have (0, 0, c, d) = (0, 0,n,n + 2) for some n G Z. Reduction by the Apollonian group 
changes (0, 0, n, n + 2) (0, 0, n, n — 2), and therefore there are only two packings containing 
a curvature zero circle, those are the packing containing (0, 0, 0, 2) and the packing containing 
(-1,0,0,1). □ 

To understand the geometric nature of the exceptional hyperbolic packings, we identify 
the real hyperbolic plane EI with the interior of the unit disk, with ideal boundary the circle 
of radius 1 centered at the origin. Recall that the curvature of a hyperbolic circle is given by 
coth r, where r is the hyperbolic radius of the circle. For circles contained in the real hyperbolic 
plane, r is positive and 

cothr = > 1 

with equality if and only if r = oo. However we will allow hyperbolic Apollonian circle packings 
which lie in the entire plane. We have a second copy of the hyperbolic plane which is the exterior 
of the closed unit disk, and we use the convention that (positively oriented) circles entirely 
contained in the exterior region have the sign of their curvature reversed, in which case coth r 
ranges from —1 to — cxd. Euclidean circles that intersect the ideal boundary (unit circle) are 
treated as hyperbolic circles of pure imaginary radius, and these cover the remaining range 
— 1 < cothr < 1. A hyperbolic circle of curvature is a hyperbolic geodesic, which corresponds 
to a Euclidean circle which intersects the ideal boundary (unit circle) at right angles, and a 
hyperbolic circle of curvature 1 is a horocycle, i.e., a circle tangent to the ideal boundary. 

The two exceptional integral hyperbolic packings are pictured in Figures 0] and El respec- 
tively. The pictures are not drawn to the same scale; in each one the dotted circle represents 
the ideal boundary of hyperbolic space, and is not a circle belonging to the packing. The 
hyperbolic packing (0, 0, 0, 2) contains exactly three circles assigned hyperbolic curvature zero, 
while the hyperbolic packing (—1,0,0,1) contains infinitely many circles assigned hyperbolic 
curvature zero. 

4. Distribution of Integer Root Quadruples 

In this section we restrict to the case k = 4m. We enumerate root quadruples v = (w, x, y, z), 
and without essential loss of generality we only consider such quadruples with L{v) = w + 
X + y + z > 0. The root quadruples are classified according to the size of their minimal 
element, say w = —n, where we suppose n > 0. In the spherical and hyperbolic cases k = ±4 
this corresponds geometrically to characterizing the distinct integral Apollonian packings that 
have an outer enclosing circle of signed curvature —n. 

Definition 4.1. Let Nroot{k; —n) count the number of integer root quadruples v = {w, x, y, z) 
to (5x)(v) = k with L{v) = w + x + y + z>0 and w = —n. 

In [HI Sec. 4] the number A'^i.oot(0, — n) of Euclidean root quadruples was interpreted as a 
class number for binary quadratic forms under a GL{2, Z) action. This interpretation was dis- 
covered starting from an exact formula for the number of (Euclidean) root quadruples having 
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Figure 4: The exceptional hyperbolic Apollonian circle packing (0,0,0,2) 




Figure 5: The exceptional hyperbolic Apollonian circle packing (-1,0,0,1) 
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smallest element — n, originally conjectured by C. L. Mallows and S. Northshield, and indepen- 
dently proved by Northshield Below we show that a similar class number interpretation 
holds generally for the case (5x)(a) = 4m. 

We write a binary quadratic form of even discriminant A = 4B^ — 4AC as 

[A, B, C] := AT"^ + 2BTU + CU^. (4.1) 

Here we follow the classical notation of Mathews jl3j . to maintain consistency with the notation 
of |S]. (The convention of Buell writes this form as [A,2B ,C], and allows odd middle 
coefficients.) A form is definite if A < 0, and it is then positive if it only represents positive 
values (and zero). Positive definite forms are those definite forms with both A and C positive, 
which follows from -4AC < 4(5^ - AC) < 0. 

The standard action on binary forms is an SL{2, Z)-action. A reduced form is one that 
satisfies 

2\B\ <A<C. (4.2) 

Every positive definite form is equivalent (under the 5'L(2, Z)-action) to at least one reduced 
form. All reduced forms are S'L(2, Z)-inequivalent except for [j4,S,A] = [A, — ii?,74] and 
[A,A,C]^[A,-A,C]. 

There is also a GL{2, Z)-action on the space of binary quadratic forms. It defines an action 
of on binary forms, sending [A, B, C] to [A, —B, C]. We define the GL{2, Z)-reduced 

forms to be the subset of 5L(2, Z)-reduced forms satisfying 

<2B < A<C. (4.3) 

Then every positive definite form is GL(2, Z)-equivalent to a unique GL(2, Z)-reduced form. 

A form is primitive if gcd(^,2i?,C) = 1 and is imprimitive otherwise. Define h^{A) to 
be the number of (primitive or imprimitive) GL{2, Z)-reduced forms of discriminant A. Let 
h{A) be the number of (primitive or imprimitive) /S'L(2, Z)-equivalence classes of forms of 
discriminant A. Since each GL(2, Z)-reduced form [y4,i?,C] is equivalent to the S'L(2,Z)- 
reduced form [A, —B, C], there holds 

/?(A) = i()i(A)+S(A)), (4.4) 

in which a(A) counts the number of (primitive or imprimitive) ambiguous classes. Here an 
ambiguous class is an S'L(2, Z)-equivalence class containing an ambiguous form, which is a 
(primitive or imprimitive) reduced form for which 2B divides A, i.e., reduced forms of shape 
[A,0,C] or [2B,B,C], cf. Mathews [13 Art. 159]. 

The usual class number h{A) is defined as the number of ( 2, Z) -equivalence classes of 
primitive forms. Similarly we let a(A) count the number of primitive ambiguous classes. If 
[A, B, C] is a primitive form of discriminant A, then for i > 1, [lA, IB, iC] is an (imprimitive) 
form of discriminant £^A. We have 

^(^) = E 4^) ' (4.5) 



£2 1 A 
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where we use the convention that h{n) = if n is not a vahd discriminant. In the same vein 
we have 



«(^) = E « (^) ■ (4-6) 

For A = —Am < one has a(A) = 2*"'"^, where t is the number of distinct prime divisors of m; 
see Mathews ^ Art. 156]. 

The following result gives a class number criterion counting root quadruples for an arbitrary 
Descartes equation Qjy = 4m. This result applies even when the associated discriminant A 
is non-negative, provided the definition of reduced form is given by (|4.2j) (resp. (|4.3|) ). though 
this no longer corresponds nicely to S'L(2, Z)-equivalence (resp. GL(2, Z)-equivalence). 

Theorem 4.1. Let n > be fixed. 

(1) The ordered root quadruples v = {w,x,y,z) with (^©(v) = 4m, for fixed m, hav- 
ing smallest element w = —n are in one-to-one correspondence with reduced integral binary 
quadratic forms [A,B,C] = AX^ + 2BXY + CY^ (primitive or imprimitive) having dis- 
criminant A = — 4n^ + 4m and a non-negative middle coefficient. If v = (w, x, y, z) with 
w = —n < X < y < z then the correspondence is given by 



[A,B,C] 



1. X " 

-n + X, -jV-n + x + y- z),-n + y 



(4.7) 



(2) When A = — 4n^ + 4m < 0, the number of root quadruples to Qx> = 4m with least 
element —n is given by the GL{2,7j)- class number 

Nrooti^m; -n) = h^{-4{n^ - m)). (4.8) 

Proof. Recall from Theorem 13.11 that an ordered integer quadruple v = (w,x,y,z) with 
w < X < y < z and L{v) > is a root quadruple if and only if 

w-\-x + y>z>0, (4.9) 

and Theorem 13.31 shows that root quadruples always satisfy < x < y < z. We now write 
w = —n to suggest that we are mainly concerned with negative integers. However, for neg- 
ative m small positive values of w are permitted, which necessarily satisfy w < -y/|m|, by 
Theorem 13. 3r 1 ) . 

The integer solutions Q-piw, x,y, z) = 4m with w = —n are in one-to-one correspondence 
with integer representations of — m by the ternary quadratic form 

QriX, Y, Z) := XY + XZ + YZ, 

where 

{X, Y, Z) := (]^{w + X + y - z),]^{w + X - y + z),]^{w - X + y + z) 
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The congruence condition w + x + y + z = (mod 2) for integer quadruples implies that 
X,Y,Z £ Z. This map is a bijection: given an integer solution to Qt{X,Y, Z) = in? — m, we 
write 

{w, X, y, z) := {-n, n + X + Y,n + X + Z,n + Y + Z) 

and an algebraic calculation shows that Qj:,{w,x,y, z) = Am. 

The root quadruple conditions (|4.9|) above translate to the inequalities 

Q <X <Y < Z. (4.10) 

Next, for any integer M, the integer solutions Qt{X^Y, Z) = M are in bijection with 
integer representations of — 4M by the determinant ternary quadratic form 

Q/^{A,B,C) := -AAC 

A solution (X, Y, Z) gives a solution Qa(^5 B,C) = —AM under 

{A,B,C) := {X + Y,2X,X + Z), 

and the inverse map is 

{X,Y,Z) := {B/2,A- B,C - B). 

Note that since —AM = B^ — AAC is even, then B must also be even, so B = B /2 must be an 
integer. 

One sees easily that the inequalities (|4.10() are equivalent to the inequalities 

<2B < A<C. (4.11) 

Finally, we recognize that the conditions 1)4. llj) give a complete set of equivalence classes of 
integral binary quadratic forms of fixed discriminant —AM under the action of GL{2,'L). 

Therefore, taking M = n'^ — m, these two steps associate to any (ordered) integer quadruple 
(— n, X, y, z) the binary quadratic form 



[A,B,C] 



1/ 

-n + X, 2 + X + y - z),-n + y 



of discriminant D := -Aip? — m). In order for the form to be positive definite, we must have 
Ti? > m. This imposes no constraint on n if m < 0, while if m > we must have n > y^|m|. 
This form is GL{2, Z)-reduced if and only if (— n, x, y, z) is a root quadruple. 

Conversely, given a form of discriminant D, we construct a Descartes quadruple 

(w, X, y, z) := (-n, n + A,n + C,n + A + C- 2B) 

which is a root quadruple if and only if [A^ B, C] is GL{2, Z)-reduced. □ 

Theorem 14. II counts the number of root quadruples N^ootiAm; —n) with L(v) > having a 
fixed smallest element w = —n m. terms of a class number, provided n? > m. This condition 
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Tl 


jv n) 




Tl 


JV 1^ n) 


i—n\ 

JV n) 


Tl 


JV n) 


JV n) 


1 


1 


+00 


11 


6 


12 


21 


6 


24 


2 


2 


2 


12 


6 


12 


22 


12 


12 


6 


2 


3 


13 


8 


12 


23 


16 


22 


4 


6 


/I 
4 


14 


6 


12 


24 


5 


O /I 

24 


5 


4 





15 


5 


19 


25 


19 


O /I 

24 


6 


2 


6 


16 


9 


16 


26 


16 


25 


7 


5 


8 


17 


12 


18 


27 


8 


24 


8 


6 


8 


18 


10 


10 


28 


10 


28 


9 


3 


11 


19 


10 


24 


29 


14 


24 


10 


8 


9 


20 


11 


20 


30 


14 


30 



Table 1: A'rootli^, — n) for small n. 

is always satisfied when m < 0. For the case m = there is one excluded value n = 
corresponding to multiples of the root quadruple (0,0, 1, 1), which falls under Theorem ESt 2). 
For positive m and for the remaining values — y^|m| < if < allowed by Theorem 13.31 the 
combinatorial correspondence with certain binary quadratic forms of discriminant A = — 4(n^ — 
m) > still applies, but the notion of "reduced form" in H4.2() is not that that normally 
used for indefinite binary quadratic forms, and A'].oot(4m, — n) is not interpretable as a class 
number. For hyperbolic Descartes quadruples (m = 1), there is one allowable value n = — 1 
not interpretable as a class number. Indeed there are infinitely many root quadruples when 
n = —1 by Theorem I3.3f 2). since {k,n) = (4m^, —m) with m = 1. 

Table^presents small values of N^^^^{—n) := A^root(— 4, —n) and N^^^{—n) := A'i.oot(4, — n). 
All hyperbolic root quadruples have —n < —1, according to Theorem I3.3r i). and all spherical 
root quadruples have n < 0, as shown in the proof of Theorem I3.3f l'). The value n = is not 
shown in the table; we have A^root(O) = given by the spherical root quadruple (0, 1, 1, 2). 

The next result uses bounds for class numbers to derive an upper bound for the number of 
root quadruples. 

Theorem 4.2. For all integers m the number Nrooti^nT-', —n) of ordered integer root quadruples 
(a, b, c, d) with Qv{^) = 4m and a = —n satisfies, uniformly for n > y^|m|, 

Nrooti^m; -n) = O (n(log n)(log log n)^) . 

The implied 0-constant is independent of m. 

Proof. Since h^{A) < h{A) we have 

A'i.oot(4m; -n) < Ji{-A{n^ - m)) = /t 

P I -4{n2— m) 

in which h[D) denotes the usual class number for primitive forms under SL(2, Z)-equivalence. 
Dirichlet's class number formula gives, for primitive negative discriminants D, that 

,(«) = !^<^L(1,X). 
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where x denotes the (primitive or imprimitive) character attached to the quadratic field 
Q(\/D), and w{D) denotes the number of roots of unity in this field, which is at most six. 
This formula is valid for class numbers of non-maximal orders as well, where the L-function 
now uses an imprimitive character. For primitive characters we have the bound 

L{l,x)=0{log\D\y, 

see Davenport 4, Chap. 14], and Ramare |16j for precise estimates. For imprimitive characters, 
we must correct by a finite Euler factor 17^25(1 — involving some subset S of primes 

dividing the discriminant. As shown in jSj Theorem 4.4], this is bounded above by ripG5(^~'~p)' 
and can lead to an extra factor of size at most 0(loglog Writing D = DqS'^ with Dq a 
fundamental discriminant, we have 

h{D) = oly^\\og\Do\Y,dll(l + - 

\ d \ S p \ ^ 

= o(^A^log|Do|(loglog5)2). (4.12) 

This gives 

MD) = o(vTD|log|Z)|(loglog|I)|)2) , 

valid uniformly for all discriminants D < —4. Here we used the estimate \ gd < S log log 5. 
The result follows on taking D = — 4(n^ — m). The O-constant can be taken independent of 
m since m? — m < 2n^ under the stated hypotheses. □ 




Remarks. (1) The upper bound of Theorem 14.21 can be strengthened when m = 0. In 
that case D = — 4n^, and one obtains the stronger upper bound O (n(loglogn)^) using (|4.12|) 
with Dq = —4 and S = n?. This improved upper bound also follows from Theorem 4.2]. 
The latter result counts primitive root quadruples only and gets a bound 0(n(log log n)). 
Theorem 14.21 counts primitive and imprimitive root quadruples (with primitivity defined by 
the auxiliary binary quadratic form (|4.7jl ). Summing over the imprimitive forms produces an 
extra factor of log log n. 

(2) In |S] a lower bound (^ log logn. ) obtained for the number of primitive Euclidean 
root quadruples, which implies 

N,oot{0,-n) = n(-^ ) . 

\ log log n J 

We are unable to obtain a lower bound for of equivalent strength for general k = 4m, because 
it appears to require obtaining good lower bounds for class numbers of general imaginary 
quadratic fields, which is a difficult problem. Using the Brauer-Siegel theorem (see |111 pg. 
328]) one can obtain for fixed m and any e > a lower bound 

7Vroot(4m,-n) > n^-' 
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valid for all n > no(m, e), in which the constant no(m, e) is not effectively computable. 

We end this section by raising some questions concerning the behavior of the total number 
of root quadruples with Qviv) = k having minimal element below a given bound. We define 
the summatory function 

T 

Sik,T):= Yl N,oot{k,-n), 

n=[^/\k\\+'^ 

Theorem 14.21 gives for nonzero k = 4m that 

T 

S{4m,T)= h^{4{n'^ + m)). (4.13) 

One expects these sums of GL(2, Z)-class numbers over a quadratic sequence to grow roughly 
like T^/-^. Can the leading term in their asymptotics be determined? The problem can be 
reduced to considering similar sums of S'L(2, Z)-class numbers using (|4.4|) to obtain 

S{Am,T) = - Y /i(-4(n2 + m)) + O ("r • 2i5ifez 



n=l^y\k\\+l 

where we made use of an upper bound for the number of ambiguous forms, as follows. If t 
denotes the number of distinct prime divisors of A < then 1)4. 6() yields 

a(A) < (i(A)2*+i < 2°(55FfeFA). 

S(4: T) 

One can also ask whether ^-^ tends to a limit as T ^ oo, and if so, compare the limiting 
value to that in (|2.7|) . One expects that there are asymptotically fewer spherical root quadruples 
than hyperbolic root quadruples. However, it is not true that A^i.oot(~4, n) < A'^rootC^, n) always. 
For example, 7Vroot(-4, 32) = 23 > 20 = iVroot(4, 32). 



5. Integers Represented by a Packing: Congruence Conditions 

In this section we restrict to the cases k = ±4 corresponding to spherical and hyperbolic 
Apollonian packings. We show there are congruence restrictions on the allowed "curvatures" 
modulo 12 of circles in integral spherical and hyperbolic Apollonian circle packings. 

Theorem 5.1. In any integral spherical Apollonian circle packing, the "curvatures" of the 
circles in the packing omit exactly three congruence classes modulo 12. 

Proof. It suffices to classify the unordered Descartes quadruples in all integral packings. We 
claim that these fall into one of two eighteen-element orbits mod 12, which are O and + 6 
(mod 12), where 

= {(2, 1,1,0), (5, 2, 1,0), (5, 4, 2,1), (6, 5, 4,1), (6, 5, 5, 4), (8, 2, 1,1), 
(8, 6, 1, 1), (8, 6, 5, 1), (9, 4, 2, 1), (9, 5, 4, 2), (9, 8, 2, 1), (10, 5, 1, 0), 
(10, 5, 5, 0), (10, 5, 5, 4), (10, 8, 5, 1), (10, 9, 5, 4), (10, 9, 8, 1), (10, 9, 8, 5)}, 
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and the notation + 6 (mod 12) means to add 6 to each coordinate of each element of O and 
reduce modulo 12. 

The theorem follows from the claim, since the curvatures in O omit the classes 3, 7 and 
11 (mod 12), while + 6 omits 1, 5 and 9 (mod 12). 

To check the claim, it is easy to verify that there are 212 possible solutions (without respect 
to order) of the spherical Descartes equation ()1.4() modulo 12. These lie in 19 orbits under the 
action of the Apollonian group. However, not all of the solutions modulo 12 come from integer 
solutions. 

If a = {w, X, y, z) is the reduction modulo 12 of a solution to the spherical Descartes 
equation, then there must exist a, h, c, d such that 

2((12a + wf + (126 + xf + (12c + y)^ + {12d + zf) = 

{12{a + b + c + d) + (w + X + y + z)f - 4 (5.1) 

If we reduce this equation modulo 24, we get 

2{w'^ + x'^ + y'^ + z"^) = {w + X + y + zf - 4 (mod 24), 

and we see that a must also be a Descartes quadruple mod 24. This eliminates 108 of the 212 
possibilities and 11 of the 18 orbits. Furthermore, if we reduce (|5.1|) modulo 48, we get the 
condition that 

2{w'^ + + + z'^) + 24(a + b + c + d){w + x + y + z) 

= {w + X + y + zf - 4 (mod 48). 

But since w + x + y + z = (mod 2), this is 

2{w^ + x^ + y^ + z^) = {w + x + y + zf -4 (mod 48). (5.2) 

Thus {w, X, y, z) must also be a quadruple modulo 48, ruling out 68 of the 104 remaining 
quadruples. Finally, a short computer search turns up integer spherical quadruples in all 36 
classes. □ 

Theorem 5.2. In any integral hyperbolic Apollonian circle packing, the "curvatures" of the 
circles in the packing omit at least three congruence classes modulo 12. 

Proof. We classify the (unordered) Descartes quadruples in an integral hyperbolic packing, 
showing that they fall in one of 81 equivalence classes modulo 12. We claim that the integer 
quadruples lie in 15 orbits under the action of the Apollonian group, given by 03,03 + 3, O3 — 
3, -03,-03 + 3,-03-3, 04,04 + 6, -04,-04 + 6, 07,07 + 3,07 + 9andOi3,Oi3 + 6, 
where 

03 = {(2, 4, 8, 8), (2, 8, 8, 8), (8, 8, 8, 10)} 
O4 = {(1, 2, 2, 11), (2, 2, 5, 11), (2, 2, 5, 7), (2, 5, 10, 11)} 
O7 = {(0, 0, 0, 10), (0, 0, 0, 2), (0, 0, 2, 4), (0, 0, 4, 6), (0, 0, 6, 8), (0, 0, 8, 10), (0, 4, 6, 8)} 
Oi3 = {(0, 0, 1, 11), (0, 0, 1, 3), (0, 0, 3, 5), (0, 0, 5, 7), (0, 0, 7, 9), (0, 0, 9, 11), (0, 1, 3, 8), 
(0, 3, 4, 5), (0, 3, 4, 9), (0, 3, 8, 9), (0, 4, 9, 11), (0, 7, 8, 9), (3, 4, 8, 9)}. 
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It is easy to check that each of these orbits omits at least three residue classes modulo 12; for 
example O3 omits the eight classes 0,1,3,5,6,7,9, and 11 (mod 12), so the theorem follows 
from the claim. 

To prove the claim, we first calculate there are 278 solutions to the hyperbolic Descartes 
equation mod 12, and then check how many of these lift to global integer solutions of the 
hyperbolic Descartes equation. Of these solutions, 184 can be eliminated by consideration 
modulo 24 and 48 as in the spherical case. However, 13 quadruples not on the above list work 
modulo 48. These 13 quadruples lie in three orbits, namely 

O3 + 6 = {(2, 2, 2, 4), (2, 2, 2, 8), (2, 2, 8, 10)}, 
O3 + 8 = {(2, 4, 10, 10), (4, 10, 10, 10), (8, 10, 10, 10)}, (5.3) 
07 + 6 = {(0, 2, 6, 10), (0, 2, 6, 6), (0, 6, 6, 10), (2, 4, 6, 6), (4, 6, 6, 6), (6, 6, 6, 8), (6, 6, 8, 10)}. 

To eliminate these quadruples, we reduce the hyperbolic Descartes equation modulo 96, which 
yields 

2((12a + wf + (126 + xf + (12c + y)^ + {12d + zf) = 

{12{a + b + c + d) + {w + x + y + z)f + 4 (mod 96). (5.4) 

After substituting {w,x,y,z) = (2,2,2,4), (2,4,10,10), or (0,2,6,10) into (^31) we get 

480^ + 48a + 48 EE (mod 96) 

where a = a + b-tc + d. This equation has no solutions as + a + 1 = 1 (mod 2) for any a. 
We have shown that one quadruple from each of the orbits (|5.3j) is impossible, thus ruling out 
the entire orbit. 

Finally, a quick search shows that the other 15 orbits all occur as reductions of integer 
hyperbolic Apollonian quadruples. □ 

In PI, it is conjectured that in any integral Euclidean Apollonian packing, all sufficiently 
large integers not ruled out by congruence conditions modulo 24 are represented. This was 
suggested by numerical data for a few specific packings. Table 121 presents integers below 100000 
missing from the spherical Apollonian packing (0,1,1,2) grouped by residue classes modulo 
24; we give data for residue classes 0,1,2,5 (mod 24) as representative. It is not completely 
clear from this data whether one should believe that only a finite number of integers will be 
missed in each class. For certain residue classes modulo 24 the exceptions seem to be rapidly 
thinning out. However for the class (mod 24) (and also 9, 18 and 21 (mod 24)) the data 
seems equivocal. 

It seems very likely that there are no congruence restrictions on integers in spherical and 
hyperbolic circle packings for any modulus prime to 6. In 8, Theorem 6.2], it is shown 
that in any integral Euclidean Apollonian packing, the "curvatures" of the circles include 
representatives of every residue class modulo m for any modulus m that is relatively prime to 
30. We expect a similar result to hold for spherical and hyperbolic packings as well. 
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n = (mod 24) 



24 


48 


144 


168 


264 


384 


504 


528 


720 


792 


864 


960 


984 


1008 


1104 


1344 


1392 


1632 


1728 


1824 


2208 


2232 


2352 


2448 


2592 


2904 


3144 


3192 


3384 


3744 


3984 


4032 


4104 


4248 


4464 


4584 


4944 


5352 


5424 


5664 


5784 


6048 


6384 


6624 


7272 


7344 


7464 


7776 


8064 


8160 


8664 


8712 


8808 


8904 


9264 


9312 


9984 


10032 


10224 


10248 


10368 


10752 


11064 


11184 


11304 


11424 


11592 


11952 


12648 


12864 


13272 


13368 


13584 


13704 


13824 


14064 


14160 


14304 


15144 


15552 


15624 


15984 


16704 


16872 


17712 


17784 


18192 


18264 


18768 


19224 


19704 


19872 


19944 


20304 


20424 


20664 


22104 


22632 


23304 


23568 


24744 


24816 


24984 


25248 


25944 


26544 


26904 


28008 


28224 


28392 


29376 


29784 


29976 


30072 


30744 


32544 


32904 


33552 


33744 


33888 


33936 


33984 


36024 


36672 


36984 


37104 


37464 


37632 


38424 


39024 


39096 


39624 


40104 


40152 


41784 


41904 


41952 


42816 


43272 


43584 


43848 


43896 


43944 


44856 


45144 


45912 


46752 


46824 


46872 


47184 


47784 


49704 


49992 


50160 


50448 


50928 


51312 


51504 


51552 


52584 


53424 


54816 


54864 


55464 


55728 


55824 


56736 


56856 


58344 


58824 


60528 


60744 


61872 


62424 


63120 


64224 


64344 


64824 


64944 


65544 


65664 


66144 


66744 


66984 


68304 


68472 


70344 


71712 


71856 


72744 


73944 


76560 


76584 


77592 


77664 


78384 


79584 


80664 


81264 


81984 


83784 


84384 


84624 


85536 


87744 


88152 


88344 


92712 


93024 


93672 


93864 


94584 


95184 


96144 


96216 


97848 


99264 


99792 


99864 


99984 


n = 1 (mod 24) 


49 


217 


529 


553 


889 


1441 


1897 


2737 


3073 


3337 


4009 


4417 


4609 


6049 


7273 


8449 


9289 


9889 


10609 


10921 


11017 


11257 


11809 


11929 


12649 


14449 


15289 


18529 


18601 


20209 


23017 


24577 


27769 


37009 


39577 


43489 


45649 


46129 


47449 


51049 


52369 


54313 


54529 


55369 


56809 


66289 


66889 


73249 


80569 


83329 


91129 


91729 


92329 
















n = 2 (mod 24) 


434 


506 


1034 


3626 


7706 


8786 


12674 


19634 


25154 


27554 


28034 


29714 


41354 


83426 


97850 












n = 5 (mod 24) 


77 


749 


869 


893 


1301 


2189 


2261 


2429 


2573 


3317 


3509 


5789 


6077 


9437 


16469 


17789 


19589 


22589 


22709 


44549 


56909 


72869 



















Table 2: Missing curvatures below 10^ in the spherical (0, 1, 1,2) packing 
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